The Dual Meissner Effect and Magnetic Displacement Currents 
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The dual Meissner effect is observed without monopoles in quenched SU(2) QCD with Landau 
gauge-fixing. Magnetic displacement currents that are time-dependent Abelian magnetic fields act 
as solenoidal currents squeezing Abelian electric fields. Monopoles are not always necessary for the 
dual Meissner effect. A mean-field calculation suggests that the dual Meissner effect through the 
mass generation of the Abelian electric field is related to a gluon condensate (A^A®) 7^ of mass 
dimension 2. 
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The understanding of the color confinement mecha- 
nism is an important problem that is yet to be solved. 
It is believed that the dual Meissner effect is the 
mechanism 0, 0. However, the cause of the dual Meiss- 
ner effect has not yet been clarified. A possible cause is 
the appearance of magnetic monopoles after projecting 
SU(3) QCD to an Abelian U(l) 2 theory by partial gauge 
fixing Q. If such monopoles condense, the dual Meissner 
effect could explain the color confinement. In fact, an 
Abelian projection adopting a special gauge called max- 
imally Abelian gauge (MA) 0, Q leads us to interesting 



results |a 13 that support importance of monopoles. 

Now, the following question arises: What happens 
when other general Abelian projections are adopted or 
when no Abelian projection is adopted? For example, 
consider an Abelian projection diagonalizing Polyakov 
loops. Monopoles exist in the continuum limit at a point 
where eigenvalues of Polyakov loops are degenerate 0. 
However, it can be easily shown that such a point moves 
only in a time-like direction. In other words, there exist 
only time-like monopoles that do not contribute to the 
string tension |8|. Let us discuss another simple case of 
the Landau gauge. Vacuum configurations are smooth 
and monopoles arising from singularities do not exist. In 
these cases, monopole condensation does not occur. We 
have to find a more general confinement mechanism that 
is realized in QCD. 

This study shows that the dual Meissner effect in an 
Abelian sense works well even when monopoles do not 
exist. Monte-Carlo simulations of quenched SU(2) QCD 
in the Landau gauge are adopted. Instead of monopoles, 
time-dependent Abelian magnetic fields regarded as mag- 
netic displacement currents squeeze the Abelian elec- 
tric fields. The dual Meissner effect implies the dual 
London equation and mass generation of the Abelian 
electric fields that may be related to the existence of 
a dimension 2 gluon condensate. The present numer- 
ical results are not perfect since the continuum limit, 
infinite- volume limit, and gauge- independence have not 
yet been examined. These discussions use the Abelian 
components based only on the assumption that Abelian 
components are dominant in infrared QCD (Abelian 
dominance 0Q,0,E|) i however, this assumption has not 
yet been clarified. Nevertheless, the authors believe that 



the results obtained here are very interesting, since they 
show for the first time that the Abelian dual Meissner 
effect works in lattice non- Abelian QCD without requir- 
ing monopoles from a singular gauge transformation[ll|. 
The gauge adopted here is the simple one and is used 
only for the purpose of obtaining smooth configurations. 
Hence, these results suggest that the Abelian dual Meiss- 
ner effect is the actual universal mechanism of color con- 
finement. Moreover, the probable relation between the 
Abelian dual Meissner effect and the dimension 2 gluon 
condensate sheds new li ght on the importance of the 
gluon condensate [3 El 111 111 111 13 El • 

We use the following improved gluonic action found by 
Iwasakipjj as a lattice SU(2) QCD: 



S = 



(3 I Co Tr(plaquette) + C\ Tr (rectangular)^ , 



which enables us to obtain better scaling behaviors of 
physical quantities. The mixing parameters are fixed as 
Co + 8ci = 1 and c\ — —0.331. 

In order to directly measure the correlations of gauge- 
variant electric and magnetic fields, we adopt the sim- 
plest gauge, called Landau gauge, which maximizes 
E s ,/i Tr [^( s ) + ^( s )]- After the gauge fixing, we at- 
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FIG. 1: Abelian Ea and non-Abelian E electric field profiles 
in Landau gauge. W(R x T = 6 x 6) is used. 



FIG. 2: Magnetic displacement currents in Landau gauge as a solenoidal current. 



tempt to measure electric and magnetic flux distribu- 
tions by evaluating correlations of Wilson loops and field 
strengths. For comparison, we also discuss the MA gauge 
where J2 S u ^' r [^( s ) cr 3^'i( s ) (J 3] i s maximized. In order 
to obtain a good signal to noise ratio, the APE smearing 
technique[20j is used when evaluating the Wilson loops 
W(R,T) = W Q + iW a a a . 

The measurements of the string tension are used to 
fix the scale when we use yfa — 440 [Me V]. We adopt a 
coupling constant /3 = 1.2 in which the lattice distance 
a(/3 = 1.2) is 0.0792(2) [fm]. This lattice distance is cho- 
sen based on the comparison between our results and 
those of the MA gauge [2ll E^ ]. The lattice size is 32 4 
and after 5000 thermalizations, we have taken 5000 ther- 
malized configurations per 100 sweeps for measurements. 

Non-Abelian electric and magnetic fields are defined 
from lxl plaquette U^ v (s) = U° v + iU^<7 a , similar to 
Ref.JH: 

E a k {s) = \{U a ik { S -k)+Ui k {s)), 

B a k {s) = l eklm (U? m (s-l-m) 

+U? m (s - 1) + U? m {s - m) + U? m (s)). 

We also define the Abelian electric (E Ai ) and magnetic 
fields {B%f) in a similar manner using the Abelian pla- 
quettes 0" {s) defined with the link variables 0"(s); 

9^( S ) = e;( s ) + e a v { s + p) - 9^ + &)- e a v { s ), (1) 

where 6^(s) is given by U^s) = exp(i6^(s)<7 a ). In the 
MA gauge, the Abelian link variable 0ff A (s) is defined by 
a phase of the diagonal part of a non- Abelian link field: 

Ufa) = ^l-Ms)\Hin9ff A (s). 

Since the off-diagonal part |c M (s)| is small@, 9ff A (s) ^ 
9^(s) in MA gauge. As a source corresponding to a static 



quark and antiquark pair, we adopt only non-Abelian 
Wilson loops in this study. 

First, we show Abelian and non-Abelian electric field 
profiles around a quark pair in the Landau gauge in Fig^ 
The profiles are studied mainly on a perpendicular plane 
at the midpoint between the two quarks. It must be 
noted that the electric fields perpendicular to the QQ 
axis are found to be negligible. It is very interesting to 
observe from Fig^that the Abelian electric field Ea z de- 
fined here is also squeezed alth oug h short-range Coulom- 
bic contributions are different [2J|. It is, thus, essential 
to ascertain the reason for the squeezing of the Abelian 
flux. 

Let us now discuss only the flux distributions of 
Abelian fields. It is numerically verified|2j| that no 
DeGrand-Toussaint monopoles|2g are present. Hence, 
the Abelian fields satisfy the simple Abelian Bianchi iden- 
tity kinematically, as demonstrated below: 

Vx£] = d 4 B%, V • B a A = 0. (2) 

In the case of MA gauge, there exist additional monopole 
current (k, k^j contributions: 

V x E MA = d A B MA + k, V • B MA = k A . (3) 

Here, E MA and B MA are defined in terms of plaquette 
variables 9jf v A (s) (mod 2n) that are constructed with 
9ff A (s). 

The Coulombic electric field arising from the static 
source is expressed in terms of the gradient of a scalar 
potential in the lowest perturbation theory. Hence, it 
contributes neither to the curl of the Abelian electric 
field nor to the Abelian magnetic field in the Abelian 
Bianchi identity, Eq.(J5J. According to the dual Meissner 
effect, the squeezing of the electric flux occurs due to the 
cancellation of Coulombic electric fields and those due 
to the solenoidal magnetic currents. In the case of the 
MA gauge, magnetic monopole currents k serve as the 
solenoidal current 0, 122I [27j ■ 

What happens in a smooth gauge like the Landau 
gauge where monopoles do not exist? It can be seen 
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FIG. 3: Curl of Abelian electric fields and magnetic displace- 
ment currents around a static quark pair in Landau gauge. 
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FIG. 4: Curl of Abelian electric fields, monopole currents, and 
magnetic displacement currents around a static quark pair in 
MA gauge. 



from Eq.J2J that only O^Ba, which is regarded as a 
magnetic displacement current, can play the role of the 
solenoidal current. It is very interesting to see Fig|3 
that demonstrates the occurrence of this phenomenon 
in Landau gauge. Note that the solenoidal current is 
in a direction that squeezes the Coulombic electric field. 
Let us also see the R dependence shown in Fig|3] The 
components of the magnetic displacement current 3\Bat 
and O^Baz do not vanish; however, they are extremely 
suppressed. In comparison, we present the case of the 
MA gauge in Fig0] Here, B^Ba^ is found to be nu- 
merically negligible as already expected from the cited 
literatures [HJ, l2ll |27fl. Instead monopole currents are 
found to circulate|2ll I22I l27j . In this case, k r is non- 
vanishing, although it is also suppressed in comparison 
with kgj,. k z is almost zero. The authors believe that 
the non- vanishing of the radial and z components of d^B 
in the Landau gauge and k r in the MA gauge is due to 



FIG. 5: Profiles of non- Abelian electric field E% after cooling. 
E a Az represents the Abelian electric field. W(R x T = 6 x 6) 
is used. 



lattice artifacts and the small size of the Wilson loop 
used here. It is interesting that the shapes of ftt-B^ in 
the Landau gauge and k^ in the MA gauge appear to be 
similar even though the strengths are different. These 
shapes have a peak at almost the same distance at ap- 
proximately 0.2[fm] and almost vanish at approximately 
0.7[fm]. 

It can be speculated whether the above consideration 
of Abelian fields is sufficient to understand the non- 
perturbative confinement problem in infrared QCD in 
terms of Abelian quantities 0, 0, HU- Let us attempt to 
verify the Abelian dominance in the Landau gauge using 
a controlled cooling 29j under which the string tension 
remains non-vanishing. It is found that the profile of the 
non- Abelian electric field E% tends to that of the Abelian 
one of E\ z in the long-range region, as shown in FigEJ 
This is consistent with the previous result [Toj using a 
different approach. 

It has been shown that the magnetic displacement cur- 
rents are important in the dual Meissner effect when 
there are no monopoles. In such a case, how can the 
origin of the dual Meissner effect without monopole con- 
densation be understood? The Abelian dual Meissner 
effect indicates the massiveness of the Abelian electric 
field as an asymptotic field: 



0. 



(4) 



This leads us to a dual London equation that is a key to 
the dual Meissner effect. Let us evaluate the curl of the 
magnetic displacement current. Using Eq.|J2J, we obtain 

Vx9A = V(V ■ E A ) -V 2 E A - 

From Eq.||lJ, we obtain the dual London equation: 

V x d 4 B A ~ (<9f - m 2 )E A - (5) 



Let us use the simple mean-field approach developed by 
Fukuda[3(j. By neglecting the gauge-fixing and Fadeev- 
Popov terms, we obtain the equation of motion D^fF^ = 



4 



and the (non-Abelian) Bianchi identity D ab *F b „ = 0. 
By applying D operator to the Bianchi identity and us- 
ing the Jacobi identity and the equations of motion, we 
obtain (D 2 p ) ab F b u = 2ge abc F b a F^ a . Note that (D 2 p ) ab = 
d 2 p S ab + ge acb {d p A c p ) + g 2 (A a p A b - S ab {A L p ) 2 ). Hence, if 
(A^Al) — 5 ab 5^v 2 ^ 0, we see asymptotically that the 
electric fields become massive {d 2 — m 2 )E% ~ with 
m 2 = 8g 2 v 2 . Now, the Abelian electric field is also mas- 
sive asymptotically (d 2 — m 2 )E\ k ~ 0. Hence, the dual 
London equation (JSJ is obtained. 

The importance of the dimension 2 gluon condensate 
has been stressed by Zakharov and his collaborators^^ 
and also by many authors in Refs. 13, 14]. Recent discus- 
sions on the value of the gluon condensate can be seen in 
Ref. 0] • Some discuss the mass generation of the gluon 
propagator [3lj . However, it is found [Isj that the gluon 



propagator has no on-mass-shell poles in the presence of 
the dimension 2 gluon condensate. This does not directly 
imply that the electric field propagator has no on-mass- 
shell pole since a gluon field A and an electric field E 
are different canonical variables The lattice study 
of gluon propagators in the Landau gauge demonstrates 
the importance of the dimension 2 condensate and its rel- 
evance to instantons|l6|. Although the operator of the 
gluon condensate is gauge variant, there have been var- 
ious discussions concerning gauge-invariant contents in 
the expectation value of the operator [l2l Il3l IbH Il7j . 

The numerical simulations of this work were performed 
using RSCC computer clusters in RIKEN. The authors 
would like to thank RIKEN for their computer facilities. 
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